Stellar Mergers Are Common

The observed Galactic rate of stellar mergers or the initiation of common envelope phases brighter than M_V=-3 (M_I=-4) is of order 0.5 (0.3)/year with 90% confidence statistical uncertainties of 0.24-1.1 (0.14-0.65) and factor of 2 systematic uncertainties. The (peak) luminosity function is roughly dN/dL L^(-1.4+/-0.3), so the rates for events more luminous than V1309 Sco (M_V=-7 mag) or V838Mon (M_V=-10 mag) are lower at r~0.1/year and 0.03/year, respectively. The peak luminosity is a steep function of progenitor mass, L M^(2-3). This very roughly parallels the scaling of luminosity with mass on the main sequence, but the transients are ~2000-4000 times more luminous at peak. Combining these, the mass function of the progenitors, dN/dM M^(-2.0+/-0.8), is consistent with the initial mass function, albeit with broad uncertainties. These observational results are also broadly consistent with the estimates of binary population synthesis models. While extragalactic variability surveys can better define the rates and properties of the high luminosity events, systematic, moderate depth (I>16 mag) surveys of the Galactic plane are needed to characterize the low luminosity events. The existing Galactic samples are only ~20% complete and Galactic surveys are (at best) reaching a typical magnitude limit of <13 mag.Comment: Submitted to MNRAS (13 pages, 6 figures, 3 tables

Constraining population synthesis models via the binary neutron star population

The observed sample of double neutron-star (NS-NS) binaries presents a challenge to population-synthesis models of compact object formation: the parameters entering into these models must be carefully chosen so as to match (i) the observed star formation rate and (ii) the formation rate of NS-NS binaries, which can be estimated from the observed sample and the selection effects related to the discoveries with radio-pulsar surveys. In this paper, we select from an extremely broad family of possible population synthesis models those few (2%) which are consistent with the observed sample of NS-NS binaries. To further sharpen the constraints the observed NS-NS population places upon our understanding of compact-object formation processes, we separate the observed NS-NS population into two channels: (i) merging NS-NS binaries, which will inspiral and merge through the action of gravitational waves within $10$ Gyr, and (ii) wide NS-NS binaries, consisting of all the rest. With the subset of astrophysically consistent models, we explore the implications for the rates at which double black hole (BH-BH), black hole-neutron star (BH-NS), and NS-NS binaries will merge through the emission of gravitational waves.Comment: (v1) Submitted to ApJ. Uses emulateapj.cls. 8 pages, 7 figures. (v2) Minor textual changes in response to referee queries. Substantial additions in appendicies, including a detailed discussion of sample multidimensional population synthesis fit

A New Template Family For The Detection Of Gravitational Waves From Comparable Mass Black Hole Binaries

In order to improve the phasing of the comparable-mass waveform as we approach the last stable orbit for a system, various re-summation methods have been used to improve the standard post-Newtonian waveforms. In this work we present a new family of templates for the detection of gravitational waves from the inspiral of two comparable-mass black hole binaries. These new adiabatic templates are based on re-expressing the derivative of the binding energy and the gravitational wave flux functions in terms of shifted Chebyshev polynomials. The Chebyshev polynomials are a useful tool in numerical methods as they display the fastest convergence of any of the orthogonal polynomials. In this case they are also particularly useful as they eliminate one of the features that plagues the post-Newtonian expansion. The Chebyshev binding energy now has information at all post-Newtonian orders, compared to the post-Newtonian templates which only have information at full integer orders. In this work, we compare both the post-Newtonian and Chebyshev templates against a fiducially exact waveform. This waveform is constructed from a hybrid method of using the test-mass results combined with the mass dependent parts of the post-Newtonian expansions for the binding energy and flux functions. Our results show that the Chebyshev templates achieve extremely high fitting factors at all PN orders and provide excellent parameter extraction. We also show that this new template family has a faster Cauchy convergence, gives a better prediction of the position of the Last Stable Orbit and in general recovers higher Signal-to-Noise ratios than the post-Newtonian templates.Comment: Final published version. Accepted for publication in Phys. Rev.

Polar kicks and the spin period - eccentricity relation in double neutron stars

We present results of a population synthesis study aimed at examining the role of spin-kick alignment in producing a correlation between the spin period of the first-born neutron star and the orbital eccentricity of observed double neutron star binaries in the Galactic disk. We find spin-kick alignment to be compatible with the observed correlation, but not to alleviate the requirements for low kick velocities suggested in previous population synthesis studies. Our results furthermore suggest low- and high-eccentricity systems may form through two distinct formation channels distinguished by the presence or absence of a stable mass transfer phase before the formation of the second neutron star. The presence of highly eccentric systems in the observed sample of double neutron stars may furthermore support the notion that neutron stars accrete matter when moving through the envelope of a giant companion.Comment: To appear in the proceedings of "40 Years of Pulsars: Millisecond Pulsars, Magnetars, and More", August 12-17, 2007, McGill University, Montreal, Canad

Double Compact Objects III: Gravitational Wave Detection Rates

The unprecedented range of second-generation gravitational-wave (GW) observatories calls for refining the predictions of potential sources and detection rates. The coalescence of double compact objects (DCOs)---i.e., neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS), and black hole-black hole (BH-BH) binary systems---is the most promising source of GWs for these detectors. We compute detection rates of coalescing DCOs in second-generation GW detectors using the latest models for their cosmological evolution, and implementing inspiral-merger-ringdown (IMR) gravitational waveform models in our signal-to-noise ratio calculations. We find that: (1) the inclusion of the merger/ringdown portion of the signal does not significantly affect rates for NS-NS and BH-NS systems, but it boosts rates by a factor $\sim 1.5$ for BH-BH systems; (2) in almost all of our models BH-BH systems yield by far the largest rates, followed by NS-NS and BH-NS systems, respectively, and (3) a majority of the detectable BH-BH systems were formed in the early Universe in low-metallicity environments. We make predictions for the distributions of detected binaries and discuss what the first GW detections will teach us about the astrophysics underlying binary formation and evolution.Comment: published in ApJ, 19 pages, 11 figure

Eccentricities of Double Neutron Star Binaries

Recent pulsar surveys have increased the number of observed double neutron stars (DNS) in our galaxy enough so that observable trends in their properties are starting to emerge. In particular, it has been noted that the majority of DNS have eccentricities less than 0.3, which are surprisingly low for binaries that survive a supernova explosion that we believe imparts a significant kick to the neutron star. To investigate this trend, we generate many different theoretical distributions of DNS eccentricities using Monte Carlo population synthesis methods. We determine which eccentricity distributions are most consistent with the observed sample of DNS binaries. In agreement with Chaurasia & Bailes (2005), assuming all double neutron stars are equally as probable to be discovered as binary pulsars, we find that highly eccentric, coalescing DNS are less likely to be observed because of their accelerated orbital evolution due to gravitational wave emission and possible early mergers. Based on our results for coalescing DNS, we also find that models with vanishingly or moderately small kicks (sigma < about 50 km/s) are inconsistent with the current observed sample of such DNS. We discuss the implications of our conclusions for DNS merger rate estimates of interest to ground-based gravitational-wave interferometers. We find that, although orbital evolution due to gravitational radiation affects the eccentricity distribution of the observed sample, the associated upwards correction factor to merger rate estimates is rather small (typically 10-40%).Comment: 9 pages, 8 figures, accepted by ApJ. Figures reduced and some content changed, references adde

Eccentric double white dwarfs as LISA sources in globular clusters

We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clusters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) mission and distances as far as the Large Magellanic Cloud, multiple harmonics from eccentric DWDs can be detected at a signal-to-noise ratio higher than 8 for at least a handful of eccentric DWDs, given their formation rate and typical lifetimes estimated from current cluster simulations. Consequently the association of eccentricity with stellar-mass LISA sources does not uniquely involve neutron stars, as is usually assumed. Due to the difficulty of detecting (eccentric) DWDs with present and planned electromagnetic observatories, LISA could provide unique dynamical identifications of these systems in globular clusters.Comment: Published in ApJ 665, L5